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Abstract: We are interested in the inverse problem of the determination of the potential 
p(x), x G SI C W 1 from the measurement of the normal derivative d v u on a suitable part T 
of the boundary of SI, where u is the solution of the wave equation duu{x,t) — Au(x,t) + 
p(x)u(x,t) = set in £1 x (0,T) and given Dirichlet boundary data. More precisely, we will 
prove local uniqueness and stability for this inverse problem and the main tool will be a global 
Carleman estimate, result also interesting by itself. 
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1 Introduction and main result 

O . 

Let n G N, T > and let SI C R™ be a bounded domain with C7 2 -boundary dSl. Let L be an 
open subset of dSl. Throughout this paper, for a functional v — v(x, t) with x G SI, t G (0, T), 
we use the following notations : 

o 

V ( ^ V ^ V \ D 2 ( ^ V \ A d ^ V 

\dxi , "" > dx n ) ' V ) x<i tj < n ' V ~[ dx V tV dt ' 

dv 

' v G M™ denotes the unit outward normal vector to <9f2 and d v v = — = Vv • 

r.^ ' du 



dv 

We consider the wave equation : 

d 2 y{x, t) - Ay{x, t) + q{x)y(x, t) = g(x, t), x G Si, t G (0, T), 
y(x, t) = h{x, t), x G dSl, t G (0, T), (1) 

y(z,0) = y°{x), d t y(x,0)=y 1 (x), x e Si. 

First of all, assuming that y° G £ 2 (0), y 1 G p G £°°(fi), /l G L 2 (<90 x (0,T)) and 

3 G L 1 (0 , T ; L 2 (SI)) are known, and assuming the compatibility condition h(x, 0) = yo(x) for 
all a; G 90, the Cauchy problem (JTJ) is well-posed and one can also prove (using a method by 
transposition, since we have non-homogeneous boundary conditions) that 

ue C{[0,T},L 2 {Sl))nC 1 {[0,T},H-\Sl)). 
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This result can be read in [10] for instance. One will also find a classical existence and regu- 
larity result when the boundary data h is equal to in Lemma Q] (see also [3]), useful in the 
inverse problem result. 

This paper treats at the same time two kinds of inverse problems which can be stated as 
follows. 

Non linear inverse Problem : Is it possible to retrieve the potential q — q(x), x € O 
from measurement of the normal derivative d u y\ roX( ^ Q T ^ where y is the solution to ([T]), Tq is a 
part large enough of the boundary <9f2 and the observation time T is also large enough ? 

We will actually give local answer to this question. If we denote by y[p] the weak solution 
of equation ([T]), assuming that p £ L°°(f2) is a given potential, we are concerned with the 
stability around p. That is to say p and y[p] are known while q is unknown and we prove the 
following local lipschitz stability result. In this direction, we will answer to two more precise 
problems. 

Uniqueness : Under geometrical conditions on Tq and T, does the equality d v y\q\ = d v y\p\ 
on To x (0, T) imply q = p on 51 ? 

Stability : Under geometrical conditions on To and T, is it possible to estimate \\q — p\\l 2 (q.) 
or better, a stronger norm of {p — q), by a suitable norm of d v y[q] — d v y[p] on T x (0, T) ? 

We will actually work on a linearized version of the inverse problem and consider the 
following wave equation : 

d?u(x, t) - Au(x, t) + q(x)u(x, t) = f(x)R(x, t), x £ O, t £ (0, T), 
u(x, t)=0, x £ dtt, t £ (0, T), (2) 

u(x,0)=0, d t u(x,0)=0, i6!!. 

Linear inverse problem : Is it possible to determine f(x), x £ ft from the knowledge 
of the normal derivative d u u\ roX ^ T ^ where R and q are given and u is the solution to ([2])? 

These questions for the wave equation have all already received positive answers since the 
uniqueness result for the linear inverse problem has been proved by M.V. Klibanov in [5] and 
Lipschitz stability results (for both linear and non-linear inverse problems) of M. Yamamoto, 
deriving from it, can be read in pTj. The proof in [14] is based on a local Carleman estimate 
for the wave operator and a compactness-uniqueness argument in order to conclude to the 
stability from the uniqueness result and an observability estimate. We aim in our document at 
giving a direct proof of a Lipschitz stability estimate from a global Carleman estimate, result 
also interesting by itself. Another gain of this new proof of the precise result given below is 
the weakened assumptions on the solution of the wave equation under study. Besides, from 
the Carleman estimate we prove in the sequel, we directly obtain that a measurement of the 
flux of the solution on a suitable part Tq of the boundary (instead of the whole boundary dfl) 
is sufficient. 

To precisely state the results we will prove in this article, we introduce, for m > 0, the set 
L^ m (Q) = {q£ L°°(G), s.t. ||g|U- (n) < to}. 
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Moreover, we also specify the geometrical assumption: 



3x ^ such that T D {x £ d£l; (x — x Q ) ■ v(x) > 0} 
T > sup |x — a?o I 



(3) 
(4) 



Theorem 1. Let m > 0, K > and r > 0. Let p belong to L^ m (Vt). Assume that the solution 
y[p] of equation ([l} is such that 



< K 



h[p]\\m(o,T-L^(Q.)) 
and assume also that the initial datum y° satisfies 

w£{\y°(x)\,xe (ft)} > r. 



(5) 



(G) 



Under the hypothesis that To satisfies the geometrical condition §3§ and T satisfies @, then 
for all q £ L^ m (il), d v y\q] — d v y[p] belongs to i? 1 (0, T; L 2 (Tq)) and there exists a constant 
C = C(m, T, K,r)>0 such that for any q e L^ m (Q), 



dy[q] dy[p] 



dv 



dv 



W(0,T;L*(r )) 



< C\\p- q\\ L 2 (n) , 



\W-p\\l'(p) < c 



dy[q] dy[p] 



di 



dv 



(7) 



(8) 



ffi(0,T;L2(r )) 



Estimate ((8j) is the Lispchtiz stability of the inverse problem while (0 gives the continuous 
dependance of the derivative of the solution with respect to the potential. 

Remark 1. The condition y[p] £ H 1 (0, T; can be guaranteed uniformly for p £ 

L2f m (f2) with more constraints on the data (y ,?/ 1 ), g, h in (TTJ). Indeed, if we assume (j/ ,^ 1 ) £ 
H*(n) x H 1 ^), g £ W hl (0,T;L 2 (fl)) and h £ H 2 (0, T\ H 2 (fl)), with the compatibility con- 
ditions h(x,0) = y (x) and d t h(x,0) — y (x) for all x £ dfl then d t y[p] solution of (fT|) 
belongs to the space C°([0, T]; n C 1 ([0, T]; L 2 (Q)), with estimates depending only on 

m and the norms of (y ,?/ 1 ), g, h in these spaces. Therefore, due to Sobolev's imbedding, 
d t y{p}£L 2 (0,T;L™(n)). " 

Remark 2. The condition on y° requires compatibility conditions for y° and h on dtt x {0}. 
In particular, \h(x, 0)| > r > 0, x £ dfl must be satisfied since \y (x)\ > r > 0, ae in ft. 

The method of proof of Theorem[T]is based on a global Carleman estimate and is very close 
to the approach of [6j , that concerns the wave equation with Neumann boundary condition and 
Dirichlet observation for the inverse problem of retrieving a potential. Actually, it also closely 
follows the approach of |14) but this work requires less regularity conditions on y. 

The use of Carleman estimate to prove uniqueness in inverse problems was introduced in [T] 
by A. L. Bukhgeim and M. V. Klibanov. Concerning inverse problems for hyperbolic equations 
with a single observation, we can refer to |12| . where the method relies on uniqueness 
results obtained by local Carleman estimates (see e.g. [4], [7]) and compactness- uniqueness 
arguments based on observability inequalities (see also [Ej). Related references [6] and [5] also 
use Carleman estimates, but rather consider the case of interior or Dirichlet boundary data 
observation. 
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2 Classical Global Carleman Estimate 



In this section, we are interested in proving a global Carleman estimate for the wave operator. 
One can find (local) Carleman estimates for regular functions with compact support in [5] , [3] 
and in |14| . 

Let us define the usual wave operator L by 

Lv = d 2 v — Av 

We consider a function v G L 2 (-T, T; H%(tt)) such that Lv G L 2 (~T, T; L 2 (fl)), and satisfying 
v(±T) = 0, d t v(±T) = on n. 

Let us now define, for x & ^, A > and j3 e (0, 1), and for all (x, i) e fi x (— T, T): 

^(a;,t) = \x-x \ 2 -fit 2 + C Q and t) = e x ^ x '^ (9) 

where Co > is chosen so that ip > 1 on O x (— T, T). We also set, for s > 0, 

to(:M) = ufotje"^*'*). 
We define the operator L p for p e (f2) by 

L p v = dfv — Av + pv 
that satisfies L p v G £ 2 (-T, T; L 2 (0)) if G L 2 (-T, T; L 2 (fl)). 

Let us first introduce the Carleman estimate we will prove by the formal calculation of 

Pw = e sv L(e- sv w). 

We have easily 

Pw = d 2 w-2s\ip(dtwd t Tp-Vw.VTP) +s 2 \ 2 ip 2 w(\d t ip\ 2 ~\VtP\ 2 ) - Aw 
-sXipw(df^j - A4i) - s\ 2 ipw(\dtij\ 2 - |VV>| 2 ) 

and if we set 

Piw = d 2 w ~Aw + s 2 \ 2 <p 2 w(\d t i(;\ 2 - |W>| 2 ) 

P 2 w = {a- l)sXtpw(d 2 <if> - Aijj) - s\ 2 ipw(\d t ijj\ 2 - |V^| 2 ) - 2s\Lp{d t wd t ^ - Vw.V?/>) (10) 
Rw = —as\ipw{d 2, 4> — Aij)) 

with a chosen later such that -pp^ < a < -^-^ (see (|16l) ). we get 

P\w + P 2 w — Pw — Rw. 

Let us now give a global (meaning "up to the boundary") Carleman inequality, following 
Imanuvilov's method [3]. One can read other versions of global Carleman estimates for hyper- 
bolic equations in [TS] and |13j . 

Theorem 2. Let us suppose that there exists xq g" £1 such that 

T D {x G dfl; (x - x ) ■ v{x) > 0} 
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Then for every m > 0, there exists Xq > 0, Sq > and a constant M — M(sq, Xq, T, m, SI, (3, Xq) 
such that for all p € L^? (f2), and for all X > Xq, s > sq: 

sX f [ e 2sv (\d t v\ 2 + \Vv\ 2 )dxdt + s 3 X 3 f [ e 2s *\v\ 2 dxdt 



< 



+ I I \Pi{e stp v)\ 2 dxdt+ f f \P 2 {e sip v)\ 2 dxdt (11) 
J-T Jo. J -t Jn 

M f [ e 2sv \L p v\ 2 dxdt + MsX [ [ <pe 2sip \d v v | 2 dadt. 
J-t Jn J-t Jr 



for all v £ H (—T,T;Hq(CI)) satisfying 

Lv e L 2 (n x (— T, T)), 
v(x, ±T) = d t v(x, ±T) = 0, Vxefl 

We do not give here extensive references about Carleman estimates for hyperbolic equations, 
but one can have a look at [3] and find references therein. 

Remark 3. Estimate pip is uniform inp whenp is in a bounded subset o/L^? m (f2). Moreover, 
this Carleman estimate is proved for any arbitrary time T . 

Proof. 

We will first prove estimate (jTTJ) with Lv in the right hand side instead of L p v. One will 
see at the end that the result hold as well for L p since p g L°°(f2 x (— T, T)). 

As we began to write, for w = e sip v we have Pw — e sip Lv and 

/ / (\Piw\ 2 + \P 2 w\ 2 ) dxdt + 2 [ [ P 1 wP 2 wdxdt 
J-t Jn J-t Jn 

= [ [ \Pw-Rw\ 2 dxdt (12) 
J-t Jn 

We will calculate and obtain a lower bound for 

T f 

/ PiwP 2 w dxdt. 

-t Jn 

The main goal of the proof will be indeed to minimize this cross-term by positive and 
dominant terms looking similar to the one of the left hand side of (|11[) and negative boundary 
terms that will be moved to the right hand side of the estimate. In the sake of clarity, we will 
devide the proof in several steps. 

Step 1. Explicit calculations 

n 

We set (Piw, -P2w) L 2(f2 X (_ T T )) = Ii,k where 1^ is the integral of the product of the 

i,fc=o 

ith-term in P\w and the fcth-term in P 2 w. We mainly use integrations by parts and the prop- 
erties of w such as w{±T) = 0, d t w(±T) = 0inf2andw = on dfl x (-T, T). 
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Integrations by part in time give easily, since d t Atp = 0, 

rT 

-t Jn 



I u = [ [ d 2 w{{a-l)sX<pw{d 2 ip - Aip)) dxdt 
J-t Jn 



r'T 



t Jn 



= (l-a)sX / / tp\d t w\ z (dti>- Aip)dxdt 
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sX 2 f [ <p\w\ 2 d 2 *p(d 2 ip - Aip) dxdt 
J-t Jn 

rT 



sA 3 [ [ <p\w\ 2 \d t ip\ 2 {d 2 ip- Aip) dxdt. 
J-t Jn 



2 

In the same manner, since dtVip — 0, one has 

In = / / d 2 w{-sX 2 yw(\d t ip\ 2 - \Vip\ 2 ))dxdt 
J-t Jn 

= s\ 2 f T f tp\d t w\ 2 (\d t ip\ 2 - \Vip\ 2 )dxdt-sX 2 f T f cp\w\ 2 \d 2 ip\ 2 dxdt 
J-t Jn J-t Jn 



-(2 + i).sA 3 f f ip\w\ 2 \d t ip\ 2 d 2 ipdxdt + ^- f f ip\w\ 2 \ViP\ 2 d 2 iPdxdt 
2 J-rJn 2 J-rJn 

-^y f_ T j n v\™\ 2 \dM\\dM 2 -m\ 2 )dxdt 



and using also integrations by part in space variable, we get 

rT 



Ji3 = / / d 2 w(-2s\<f(d t wdtip- Vw.Vip))dxdt 
J-t Jn 

= sX [ [ ^\d t w\ 2 d 2 ipdxdt + sX 2 [ [ tp\d t w\ 2 \d t ip\ 2 dxdt 
J-t Jn J-t Jn 

+sX [ [ ip\d t w\ 2 Aip dxdt + sX 2 [ [ tp\d t w\ 2 \Vip\ 2 dxdt 
J-t Jn J-t Jn 



\2 



-2sX I / <pd t w d t ipVw.Vip dxdt. 
-t Jn 



We compute in the same way 

rT 
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[ [ -Aw((a - l)sXifw(d 2 tp - Aip)) dxdt 
J-t Jn 

-(l-a)sxf [ tp\Vw\ 2 (d 2 tp- Ai>)dxdt 
J-t Jn 

+ (l^_o) sA 2 f f ^wfAipidfy- Aip) dxdt 
2 J-t Jn 

+ ( 1 ~ a h x 3 f [ v \w\ 2 \ViP\ 2 {d 2 ^- A^)dxdt 
2 J —T J n 



and 

fT 



J 22 = / / -Aw(-s\ 2 cpw(\d t ip\ 2 - \Vijj\ 2 ))dxdt 
J-t Jn 

= -s\ 2 f [ ^\\7w\ 2 {\d t ip\ 2 - \\7^\ 2 )dxdt 
J-t Jn 

,2 pT [ 

- S Ar \ \ v\w\ 2 A(\V^\ 2 )dxdt 
2 J-t Jn 

+^ I* I v \w\ 2 M{\dM 2 -\^^\ 2 )dxdt 

+ S ~^~ f T l/\ w \ 2 \^\ 2 (\d^\ 2 ~ \^\ 2 )dxdt 

~sX 3 [ [ ip\w\ 2 Vip-V(\Vip\ 2 )dxdt 
J-t Jn 

Using the fact that w|sox(-t.t) = 0, we have, on d£l x (— T,T), Vw = {d u w)v that give 
|u>| 2 = |<9i,u;| 2 . Therefore, we obtain 

I23 = [ [ -Aw(-2sX(fi(d t wd t ^ -Vw -Vip))dxdt 
J-t Jn 

= sX [ [ tp\\7w\ 2 (d 2 tp - Atp)dxdt + 2sX 2 [ [ <p\Vi/> ■ Vw| 2 dxdt 
J-t Jn J-t Jn 

-2sX 2 [ [ ipdtwdtipVw -\7tpdxdt + sX 2 [ [ tp\\7 w\ 2 (\d t ip\ 2 - |V?/>| 2 ) dxdt 
J-rJn J-rJn 

-sX [ [ (p\d u w\ 2 \7ip ■ vdadt + 2sX [ [ ipD 2 ip\\7w\ 2 dxdt 
J -t J en J-rJn 

since D 2 ip is symmetric. 
One easily write 

J31 = / / s 2 X 2 ip 2 w(\d t ip\ 2 - \Wtp\ 2 ){{a-l)sXipw{d 2 ^~ A^)) dxdt 
J-t Jn 

= (a-l)s 3 X 3 [ [ tp 3 \w\ 2 (dfy- AiP) (\d t ip\ 2 -\Vip\ 2 ) dxdt 
J-t Jn 



and 



ha = f T [ s 2 x 2 v 2 w{\d t ^\ 2 - |v^| 2 )(- S A 2 ^(|a t ^| 2 - |vv>| 2 )) dxdt 

J-t Jn 

-* 3 A 4 , 

' -t Jn 



■■ f T [ v 3 \w\ 2 {\d t ^\ 2 -\w^\ 2 ) 2 dxdt. 

J-t Jn 
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Finally, some integrations by part enable to obtain, since \7(\d t ip\ 2 ) = 0, 

J33 = / / s 2 \ 2 tp 2 w{\d t ip\ 2 -\Vip\ 2 ){-2sX<p(d t wd t ip-Vw.Vip))dxdt 
J -t Jn 

= s 3 X 3 f T f tp 3 \w\ 2 {d^-A^)(\d t ^\ 2 -\V^\ 2 )dxdt 
J-t Jn 



+s 3 X 3 f f (f 3 \w 2 \(2d 2 ^\d t ij\ 2 + W^-W(\W^\ 2 ))dxdt 
J-t Jn 

+3s 3 X i f T f <p 3 \w\ 2 (\d t ^\ 2 - \Vip\ 2 ) 2 dxdt. 
J-t Jn 



Gathering all the terms that have been computed, we get 

i-T 

P1V0P2W dxdt 



II 

J-t Jn 

= 2s\ [ [ (p\d t w\ 2 d 2 ipdxdt-as\ [ [ ip\d t w\ 2 (d 2 ip - kip) dxdt 
J-t Jn J-t Jn 

+ 2sX 2 ( [ ip (\d t w\ 2 \d t ^\ 2 -2d t wd t ipS7w • VV> + |V</> • Vw| 2 ) dxdt 
J-t Jn 

+ 2sX [ [ LpD 2 ip\\7w\ 2 dxdt + as\ [ [ ip\ Vw\ 2 (d 2 4> - Aip) dxdt 
J-t Jn J-t Jn 

— s\ / / tp |(9;yw| 2 Vip ■ v(x) dadt 
J-t Jan 

2s 3 X i f [ <p 3 \w\ 2 (\d t ^\ 2 ~ \V^j\ 2 ) 2 dxdt 
J-t Jn 



(13) 



+ 2s 3 X 4 r 

J-t Jn 

+ s 3 X 3 [ [ tp 3 \w 2 \{2d?iP\d t ip\ 2 + \7->p -\7(\\7ip\ 2 ))dxdt 
J-t Jn 

+ as 3 X 3 f [ <p 3 \w\ 2 (d 2 yj- A^){\d t ij\ 2 -\\7ip\ 2 )dxdt 
J-t Jn 

+ X 1 

where X\ satisfies, using the regularity of ip and the fact that > 1 implies A < e A ^ = tp, 

\X x \<MsX 3 ( ( tp 3 \w\ 2 dxdt. (14) 
J-t Jn 

Here and in the sequel, M > corresponds to a generic constant depending at least on Q and 
T but independant of s and A. 

Step 2. Bounding each term from below 

On the one hand, one can notice that 

2sX 2 [ ftp (\d t w\ 2 \d t ip\ 2 - 2d t wd t ^w • W + |W> ■ Vw| 2 ) dxdt 
J-t Jn 

= 2sX 2 [ ftp {dtwdtip - Vw • \7ipf dxdt > 0. (15) 
J-t Jn 



Besides, considering the terms in sA that must now give the dominant terms in |9 t w| 2 and 
|Vui| 2 and thus have to be strictly positive, one can guess that we need 

2d 2 tp - a(d 2 ip - AV>) > and 2D 2 ^ + a{d 2 ^ - Aip) > 0. 

This will constrain the value of the unspecified constant a > in the definition (JTUJ) of P2W. 
We get, by explicit computations, that (3 £ (0, 1), and a must satisfy 

W <«<^- (16) 



P + n P + n 

As a direct consequence, we can write 

2?i2„,,j„j, „,„\ / / ,jo „.,\2{a2„ 



2sA / / <p\d t w\*dfipdxdt-as\ / / vj|<9 t w| (<9 t i/> - cfecft 

-T ./n 7-T JO 

+ 2sX [ I ipD 2 ^\Vw\ 2 dxdt + asX [ [ <p\Vw\ 2 (d?i/; - Aip) dxdt 

J —T JQ J-tJq 



>MsX [ [ (p\d t w\ 2 dxdt + Ms\ [ [ <p\Vw\ 2 dxdt. (17) 

J-T Jq J-T Jq 



i-TJQ J-TJQ 

On the other hand, we can observe that: 

2s 3 A 4 [ T [ v 3 \w\ 2 {\d t *P\ 2 ~ \V^\ 2 ) 2 dxdt 
J-t Jn 

+ s 3 X 3 [ f ip 3 \w 2 \{2d 2 ip\d t ^\ 2 + -V{\V^\ 2 ))dxdt 

J-T JQ 

+ as 3 A 3 f [ cp 3 \w\ 2 (d 2 ip- A^>)(|9 t ^| 2 - \Vip\ 2 ) dxdt 

J-T Jq 

r T 

„3\3 / / ,„3|„„2i 



= s 3 A 3 / / Lp A \w z \F x {4>)dxdt 
J-t Jq 

where 

F x {ct>) = 2\{\d t ^\ 2 - |W| 2 ) 2 + {2d 2 4>\d t 4>\ 2 + • V(|W| 2 )) + a(dfy - A^j)(\d t ^\ 2 - |W| 2 ) 
= 32A(/? 2 i 2 -\x- x Q \ 2 ) 2 - 16(/3 3 i 2 - \x - x \ 2 ) - 8a((3 + n)(/3 2 t 2 - \x - x \ 2 ) 
= 32X(f3 2 t 2 -\x- x \ 2 ) 2 - 8(a(/3 + n) + 2fi){p 2 t 2 -\x- x \ 2 ) + 16(1 - (3)\x - x \ 2 . 

Since xq £ £7, we have 16(1 — f3)\x — xq\ 2 > c* > 0. Therefore, we are considering a polynome 
P(X) > 32XX 2 - 8 (a(fi + n) + 2(3) X + c* and taking A > large enough, the minimum of P 
will be strictly positive. Consequently, 

s 3 A 3 / ( Lp 3 \w 2 \F x {4>)dxdt>Ms 3 X 3 [ [ ip 3 \w\ 2 dxdt. (18) 
J-tJq J-tJq 



Thus, plugging ((T3J), <[T7J) and ATS} in ([13]), we obtain 

P1WP2W dxdt + 2s\ / / ip Iduw] (x — xq) ■ v(x) dadt — X\ 



T JQ J-TJdQ 

>Msx[ [ ip(\d t w\ 2 + \\7w\ 2 ) dxdt + Ms 3 X 3 [ [ ip 3 \w\ 2 dxdt. 
J-tJq J-tJq 
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Since we also have 

pT p t-T r pT 



[ [ \Pw - Rw\ 2 dxdt < 2 / / \Pw\ 2 dxdt + 2 / / \Rw\ 2 dxdt 

<M ( [ \Pw\ 2 dxdt + Ms 2 A 2 [ [ ip 2 \w\ 2 dxdt, 
J-t Jn J-t Jn 



using ([T2]) and ([13]). we get 



sX f f (\d t w\ 2 + \Vw\ 2 )<pdxdt + s 3 X 3 [ f \w\ 2 cp 3 dxdt 

+ f f (\P lW \ 2 + \P 2 w\ 2 ) dxdt 
J-t Jn 

<M [ [ \Pw\ 2 dxdt + MsX [ [ (p\d„w\ 2 {x - x ) ■ v(x)d<rdt 
J-t Jo. J-t J an 

+ MsX 3 [ [ ip 3 \w\ 2 dxdt + Ms 2 X 2 [ [ ip 2 \w\ 2 dxdt. 
J-t Jn J-t Jn 

We take now Sq large enough so that the terms of the last line (coming from X\ and |i?u>| 2 ) 
are absorbed by the dominant term in s 3 A 3 |u;| 2 9s 3 as soon as s > sq. Using also the condition 
© on r , we finally obtain for some positive constant M = M(sq, Xq, m, f2, /3, Xq), 

sX [ [ ip{\d t w\ 2 + \X?w\ 2 )dxdt + s 3 X 3 [ [ ip 3 \w\ 2 dxdt 
J-rJn J-rJn 

{ \Piw\ 2 dxdt+ I I \P 2 w\ 2 dxdt (19) 
rJn J-rJn 

<M [ [ \Pw\ 2 dxdt + MsX [ [ tp\d u w\ 2 dcrdt 
J-t Jn J-t Jt 

Vs > s , VA > A . 

Step 3. Return to the variable v 

Using that w = ve sip gives for all x e il and t e (-T, T) 



e 2s *\d t v\ 2 <2\d t w\ 2 +2s 2 \d t cp\ 2 \w\ 2 , 
e 2 ^|Vw| 2 < 2|V«f + 2s 2 |V^| 2 |u;| 2 , 
e 2sv \d v v\ 2 = \dM 2 on dn 



and that by construction Pw = e sip Lv, we can go back to the variable v in (|19p and obtain 
that there exists some positive constant M = M(so, Ao, m, J7, /3, xq) such that for all s > sq 
and A > Ao, 

sA ( T f e 2sv {\d t v\ 2 + \\7v\ 2 )dxdt + s 3 X 3 f f e 2sip \v\ 2 dxdt 
J-rJn J-rJn 

f \P 1 {e sv v)\ 2 dxdt+ f f \P 2 {e stp v)\ 2 dxdt (20) 
-rJn J-rJn 

<M [ [ e 2stp \Lv\ 2 dxdt + MsX [ [ ipe 2sv \d u v\ 2 dodt. 
J-rJn J-tJth 
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It concludes the proof of a Carleman estimate for the operator L = d 2 — A. 



Step 4. Wave operator with potential 

The Carleman estimate (|TT|) for the operator L p = d 2 — A + p with p £ L< L(0) is a direct 
consequence of (|20|) noticing that on SI x (— T, T), 

M 2 < 2| V| 2 + 2||p||i- m(n) M 2 < 2| V| 2 + 2toH 2 . 
Indeed, choosing so ( or -^o) large enough, one can absorb the term 

2Mm ( I e 2sip \v\ 2 dxdt 



T JO. 



in the left hand side of (|2"0)) and obtain (fTTj) with slightly different constants. This ends the 
proof of Theorem [21 

3 Stability of the inverse problem 

Before giving the proof of Theorem [TJ we will begin this section by a stability theorem for the 
linear inverse source problem stated after equation ^ . The following answer is obtained using 
the global Carleman estimate given above. 

Theorem 3. Let m > 0, K > and r > 0. Let q belong to L< m (fi). Assume that f € L 2 (Q) 
and R<E iJ 1 (0,T; J L oo (^)) with 



and 



l-R||ff 1 (o,T ; L°°(n)) < K 



inf \R(x,0)\ > r. (21) 



Under the hypothesis that Tq and T satisfy the geometrical conditions 

3xq $jL £l such that Tq D {x £ dfl; [x — x$) ■ v{x) > 0} and T > sup \x — x$\ 

xgfi 

if u[f] is the solution of equation 

d 2 u(x, t) - Au(x, t) + q(x)u(x, t) = f(x)R(x, t), x £ Q, t G (0, T) 
u(x,t) = 0, x £ dCl, t e (0,T) 

u(x, 0) = 0, d t u(x, 0) = 0, x £ fl, 

then there exists a constant C = C(T,Q,Tq, K,r) > such that for all f £ L 2 (Q) : 

du[f] 



c- l \\f\\ L * m < 



<C||/|U= ( o)- (22) 

ffi(0,T;L2(r )) 



Proof. We will apply the Carleman estimate given by (|20p to w = xdty where x is a cutoff 
function to be detailed later. We divide the proof in several steps. 

Step 1. Let us first work on the equation satisfied by z = dtu: 

d 2 z(x, t) - Az(x, t) + q(x)z(x, t) = f(x)d t R(x, t), x £ fl, t £ (0, T) 
z(x,t) = 0, x£dn, t£ (0,T) (23) 

z(x, 0) = 0, d t z{x, 0) = f(x)R(x, 0), x £ Q. 

We gather in the following lemma the classical energy and trace estimates we will need in the 
sequel. 
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Lemma 1. Assume that p G L% m (Sl), g £ L x (0, T; L 2 (S7)), u° G ff^(O) and u 1 G L 2 (17) and 
consider the classical wave equation 



d 2 v(x, t) - Av(x, t) + q(x)v(x, t) = g(x, t), x e 17, t e (0, T) 



v(x, t) = 0, 

v(x, 0) = vq(x), d t v(x,0) =vx(x)), 



x e <917, t G (0, T) 
£ G 17. 



(24) 



TTie Cauchy problem is well-posed and equation (|24D admits a unique weak solution 

v g c([o,T],ff 1 (n)) n c^M, 

and i/iere exists a constant C = C(Cl,T, m) > such that for all t G (0, T), £/ie energy 
E v {t) = \\dtv(t)\\^2(Q-) + l|Vu(t)||^ 2 (f2) °/ ^ e s y s t em satisfi 



E v (t) < C (\W\\ 2 H i {n) + ||wi||i2 ( n) + ll5llii( ,T,i 2 (f2))) ■ 
Moreover, the normal derivative d v v belongs to L 2 (0, T; L 2 (<917)) and satisfies 

< C (ll w o||^i(o) + ll u i|lz,2(n) + llffllz,i(o,T,£ 3 (n))) • 



dv 
dv 



(25) 



(26) 



L 2 (0,T;L 2 (dQ)) 



This result is very classical and we refer to [10] (Chapter 3) for the proof of the existence 
and uniqueness of solution to equation ([24)) . Estimate (|25|) can formally be deduced from 
multiplication of (|24|) by U{, and the integration of this equality on (0, T) x 17, using some 
integrations by parts. Concerning estimate (|2"6")h we refer to [9] (Chapter 1). This estimate is 
a hidden regularity result which can be obtained by multipliers technique. 

We can apply this lemma to equation ([23} since / G L 2 (17) and R G i/ 1 (0, T; L°°(17)). We 
denote the energy of the system by 

= \W)\\]il{n) + \\9tz{t)\\l Hn) 

and we get z G C([0, T],Hq(Q,)) n C 1 ([0, T], L 2 (17)) with, for all i G (0,T): 

E z (t) < CE z (0) + C\\fd t R\\l HO)T . L2(n)) 

< c\\f\\ 2 L 

and G L 2 (0,T;L 2 (<917)) with 



(27) 



<9z 



L 2 (0,T;L 2 (0n)) 



< Cll/I&cn) (||i*(Q)|li~ (Q ) 



\R\ 



h 1 (o,t-,l°°(si)) 



This last estimate gives at the end d v u G i? 1 (0, T; i 2 (9il)) and proves the right hand side of 
the two sided estimate (|22[) . It thus gives a meaning to the measurement of the flux of the 
solution u[f] we make in our inverse problem. 

Step 2. Let us now extend the problem (|23|) on (— T,T) in an even way, setting z(x,t) = 
z(x,—t) for all (x,t) e 17 x (— T, 0). We also extend dtR on an even way and keep the same 
notations for the new problem. Therefore, we have 



z G C([-T,T\,H%(fl)) nC 1 ([-T,T],L a (fi)) and d t R G L 2 (—T, T; L°°(17)). 
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The goal we have in mind is to apply the Carleman estimate to a solution of the wave 
equation in order to be able to estimate a weighted norm of / by a weighted norm of the 
measurement d u u. The solution z could be a good candidate but is not equal to zero at time 
±T. This motivates the use of a cutoff function in time, that is calibrated to work well with 
the Carleman weight function ip. 

We recall that for j3 £ (0, 1), ip and (p are defined by (j9]). Since T satisfies (j4|, one can choose 
/3 such that j3T 2 > sup xe0 \x-x \ 2 so that for all x £ ft, tp(x,±T) = \x-x \ 2 - /3T 2 + C < C . 
We then choose r] > such that 

iP(x, t) < C , V(x, i)e!lx [-T, -T + V ]U[T- r,, T]. 

Besides, one observes that 

ip(x,0) = \x - x \ 2 + C > C , VxeO. 

Therefore, for all (x, t) € Cl X [-T, -T + rj] U [T - -q, T], 

tp(x, t) < e XCo < e A ^o) = ^ o) (28) 

but we also have, for all (x,t) e fl x [— T, T], 

¥>(a:,t) <¥>(!, 0). (29) 

Using the parameter rj introduced here, we define the cut-off function \ € C°°(K; [0,1]) 
such that 

' X (±T) = x'(±T) = 

X (t) = l Vte[-T + r,,T-rj\ [M) 

and we set v — x z that satisfies the following equation: 

d 2 t v -Av + qv = xfdtR + x"z + 2 X 'd t z, fl x (0, T) 

v(x,t)=0, x e dfl, t e (0, T) (31) 

w(x,0)=0, a t u(a;,0) = f(x)R(x,0), x e fl. 

Henceforth, M > will correspond to a generic constant depending on sq, \q,T,Q, Tq, f3, 
X, t, K , and rj but independent of s > sq. 

Step 3. We use now the same notations as in the proof of the Carleman estimate. For v 
solution of (l3TT) . we set 

w = e sv v , Piw = d 2 w -Aw + s 2 \ 2 l p 2 w{\d t tlj\ 2 - |V?M 2 ), 

and we then have for all x £ £1, w(x,±T) = dtw(x,±T) = and w(x,0) = and for all 
(x, t) £ dQ x (-T, T), w(x, t) = 0. 

Inspired from an idea of O. Yu. Imanuvilov and M. Yamamoto [B], we consider the integral 

o , 

/ P\w{x,t) dtw(x,t) dxdt. 
-t Jn 
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On the one hand, using the properties of w, we can make the following calculation: 



-T Jn 




2 



P\w dtw dxdt 

(d 2 w - Aw + s 2 X 2 ip 2 w(\d t ip\ 2 - IVV'I 2 )) dtwdxdt 
\d t w(0)\ 2 dx - ^ [ f \w\ 2 d t ( V 2 (\d t ^\ 2 - |wi 2 )) dxdt 



t Jn 



~ / \d t w(0)\ 2 dx~2s 2 X 2 f I \w\ 2 d t (<p 2 (f3 2 t 2 -\x-x \ 2 )) dxdt 



-t Jn 

o 



1 j e 2 ^(°)|/J2(0)| 2 cfe + 8s 2 A 3 / / \w\ 2 tp 2 /3t(p 2 t 2 - \x - x \ 2 ) dxdt 



n j-tjq, 
o 

2\2 / / I |2 ,2o2j 



- 4s 2 X 2 / / \w\'y 2 2 t. 
J-t Jn 

On the other hand, from this equality and a Cauchy-Schwarz estimate, 

f e 2sip( - 0) \fR(0)\ 2 dx = 2 f [ P lW dtwdxdt + 8s 2 \ 2 [ [ w 2 ip 2 f3 2 t dxdt 
Jn J-t Jn J-t Jn 

-16s 2 A 3 f f \w\ 2 tp 2 j3t{(3 2 t 2 ~ \x - x \ 2 ) dxdt 
J-t Jn 

< 2 (^J J \P x w\ 2 dxdt) (^J J \d t w\ 2 dxdt^j 

+Ms 2 X 3 [ [ \w\ 2 ip 2 \j3t{(3 2 t 2 -\x-x \ 2 )\ dxdt. (32) 



t Jn 



Using now the 'intermediate" Carleman estimate (|19l) for a fixed A > Ao not mentioned any- 
more, choosing s large enough to absorb the last term in the right hand side (|3"!?|) . we obtain 



e 

n 



< -j= (^J J \P lW \ 2 dxdtj (s J J \d t w\ 2 dxdtj +Ms 2 X 3 J J \w\ 2 tp 2 dxdt 



^ m( r r \ Pw ]2 dxdt + s f ( v \Q vW ^dadt 

V s \J-rJn J-tJth j 



< f [ e 2sv \Lv\ 2 dxdt + s f [ e 2sv \d v v\ 2 dadt) (:«) 

yJ-T Jn J-t Jr j 



the weight function <p being bounded from above and below. 

From equation (|3Tj) and the properties (|30| of the cut-off function \ and that of the weight 
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ip given in ([2"5| and (|^|) , one gets that 
' e 2sip \Lv\ 2 dxdt 



-T Jn 



< Ml e 2s *\ X fd t R\ 2 dxdt + M e 2sv (\ X 'd t z\ 2 + \ X " z\ 2 ) dxdt 

J-t Jn J-t Jn 

?T r ( r -T+rj 



< m[ f e 2sip \f\ 2 \d t R\ 2 dxdt + M ( f "+ f ) f e 2sv (\d t z\ 2 + \z\ 2 ) dxdt 

J-T JO. \J-T JT-rtJJn 

< M[ f e 2s '^\f\ 2 \d t R\ 2 dxdt + All f +[ ) [ e 2seXC ° (\d t z\ 2 + \z\ 2 ) dxdt 

J-TJn \J-T JT-nl Jn 



< M\\d t R\\ LH _ TtT . iL - m [ e 2s ^\f\ 2 dx 

Jn 



in 

, , - T +v 

„2se AC 



+ Me 2se ^° I + \ (\d t z\ 2 + \z\ 2 ) dxdt. 

yJ-T Jt-t) J Jn 

Using the energy estimate given in (|27|) . and again the property (|2"8)) of the weight <p, one gets 

/ [ e 2s *\Lv\ 2 dxdt < M [ e 2s ^\f\ 2 dx + Me 2seXC ° ( [ +[ I E z (t)dt 
J-t Jn Jn \J-t Jt-v J 

< M [ e 2s ^\f\ 2 dx 

Jn 

+M V (\\R(0)\\U (n) + \\R\\ 2 HH0 , T ., L eo m ) e 2se " C ° J \f\ 2 dx 

< M [ e 2s ^\f\ 2 dx. 

Jn 

Gathering this last estimate with ([33| , we have proved 

f <pvW\f R (p)\* fa < f e 2s ^\f\ 2 dx + M^- 3 f [ e 2sv \0 u v\ 2 dadt. 
Jn v s Jn J-t Jr 

Therefore, the assumption (f2"TT) made on R allow to obtain 

/ e 2s ^\f\ 2 dx< M= I e 2s ^\f\ 2 dx + MV~s C [ e 2s ^\d„v\ 2 dadt, 
Jn V s Jn J-t Jr 

and the choice of s large enough gives 

f e 2 SV (0)^2 dx < My/s f I e 2sv \d„v\ 2 dadt 
Jn J-tJtq 

< Myfs f f e 2s,fi \d„(dtu)\ 2 dadt. 

J — T J Tn 



The proof of Theorem [3] is then complete. 
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